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SEMI-ANNUAL REPORT 


This progress report for grant NAG 5-145/Supplement 2 covers the period 
January through June 1984. Most of the research discussed here concerns the 
pole tide, the oceanic response to the Chandler wobble. The reader should 
consult earlier progress reports for relevant background material. 


During the beginning of the reporting period my efforts were devoted to 
completing the major project of Supplement 2, a re-examination of Wunsch's 
(1974) North Sea pole tide theory. Wunsch had proposed that observed pole 
tide enhancements there (i.e. larger amplitudes than a static tide would 
possess) resulted from bottom friction, with drag coefficient R, in 
combination with the depth h of the North Sea decreasing southward. However, 
his analysis was marred by errors (Wunsch, 1975) and improper approximations. 
As described in the previous status report, our re-analysis yielded one case 
where the equations appeared capable of satisfying all necessary boundary 
conditions: this was a situation in which the North Sea depth was linear in 
northward distance y, i.e. h = h Q o<y, and where the drag was inversely 
proportional to h(y); Wunsch had modeled the depth as proportional to l+«y and 
treated R as constant. We found the tide's stream function Y to be given by 
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where x is eastward distance. 
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and f is the Coriolis parameter. This stream function satisfies 3 of Wunsch's 

boundary conditions, Y = 0 (i.e. no flow) at the North Sea boundaries x=0, 

x“b, y=0. The unknown coefficients A would be determined from the fourth 

n 

boundary condition, at the open northern edge of the Sea. Wunsch specified 
this condition, involving the departure T' of the pole tide height from the 
static value, as 


T' = B Q exp(ip 0 x) at y=a , 
where p Q and Bq are constants. 


At the end of the previous reporting period, I had derived a quite 
complicated expression for T’(x,y) based on a rigorous general relation 
between T' andY; determination of the A n from application of the boundary 
condition was sure to be troublesome. At the beginning of the present 
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reporting period, I discovered a "shortcut” method of determining T’(x,y); the 
result was much simpler, namely 
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Unfortunately, determination of the A revealed fundamental difficulties 

n 

with the theory. In typical Fourier analysis problems, the sine and cosine 
coefficients are independent unknowns, determined separately by sine and 
cosine orthogonality. Our solution for T*, however, is in terms of sine and 
cosine coefficients which are related (and the relation is different for each 
n!); sine and cosine orthogonality will thus yield different values for An. 

We found 
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using cosine orthogonality, or 
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using sine orthogonality, where B f - Keal{B Q ?, B i - Imag {b« ] , and 
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These expressions are equal if 
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This may explain why Vunsch (1974) chose a complex boundary condition at y=a, 
even though the rest of his equations were strictly real-valued: the 
non-uniqueness is apparently eliminated if the phase (arg) of Bg is pre-set. 
Unfortunately, this approach is inconsistent because arg^Bg} must be a 
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function of n — i.«., for each n a different phase mist be selected for Bq. 

Such difficulties would not be eliminated if a boundary condition of the 

fora 


T' - B Q q(x) at y-a , 

for any function q(x), were specified, because after q(x) were expanded in a 

Fourier series the same ambiquity would appear — unless the Fourier expansion 

of q(x) happened coincidentally to possess the arg * - i* - <X n type of phases. 

Since the governing equation for 4* is an elliptical type of differential 

equation, its solution requires four boundary conditions [see Carrier & 

Pearson, 1978], so we cannot simply omit the troublesome boundary condition at 

y-a (we'd still have to determine the A ). We must instead conclude that 

n 

either Wunsch's theory is ill-posed, e.g. the original momentum equations mist 
Include other terms such as terms or diffusive friction or or else 

additional theory external to the North Sea is required in order to delineate 
the proper y-a boundary condition (this may be quite similar to Stoke' s 
Paradox! — see Tritton, 1977, pp. 92-93). 


Computer programs were written to calculate the solutions for and thus 
T' for 5 situations: I) A fl determined from cosine orthogonality, with B^ - 0 

artifically; II) A n determined as above with B^/B r = -tan( XN- ctf^) illegally; 

III) A determined from cosine orthogonality when the boundary condition is 

n 

T' - B r cos(p Q x) at y-a ; 

IV) same as III but with 


B* = all* */<■ )(*'-») { - /f } 
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where is .the colatitude of the North Sea at y-a and all other symbols are 
from Dickman (1983); and V) same as IV but with larger by a factor of VT. 
In the third situation the A ft are given by 
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The fourth and fifth situations correspond to forcing the North Sea pole tide 
by relative motion between the wobbling solid earth and wobbling oceanic body, 
as in Dickman (1983); this was the goal of the project! In IV time-averaged 
wobbles are used, while in V maximum wobble amplitudes are considered. 

Results are as follows: 
v -4 

I) Aj - -2.0x10 , ... , and -0.29 cm£T' 13.1 cm along the southern coast 

of the North Sea; 

II) A. ■ 5.2x10 ... , and -0.07 cm < -T' < 3.7 cm along that coast; 
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III) Aj “ -3.1x10 , ... , and -0.41 cm T* £ 20.7 cm along that coast; 

IV) Aj ■ 4.7x10 ... , and -0.06 cm <_ -T ' < 3.2 cm along that coast 

if ocean-solid earth coupling is 99. 9% effective (M'-M * 0.001 x M); and 

V) Aj * 6.6x10 **, ... , and -0.09 cm £ -T* £ 4.5 cm under the same conditions 
as in IV. We can conclude, firstly, that (because I-III results are so 
different) the ill-posed nature of the theory is a serious deficiency; and, 

<t±r+k 

secondly, that forcing by ocean-solid^wobbling at the edge of the Sea may 
Indeed explain North Sea pole tide observations (the enhancements T' in IV 
have precisely the correct magnitude range!). 

Dissipation of energy by North Sea pole tide currents was also computed; 
preliminary results indicate that such dissipation may explain a significant 
fraction of Chandler wobble energy loss. 


I was scheduled to devote four months to substantiating one aspect of 
earlier NASA-grant work. Dickman (1983) had found an explanation for the 
Markowitz wobble, that it is a natural wobble of the coupled rotating 

ocean-solid earth system, subject to the assumption that ocean - solid earth 
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coupling is primarily non-dissipative and of magnitude M0 N-m. The 

four-month project was to determine, using fluid dynamics, the topographic 

drag exerted on a continent by pole tide - type currents. This was, however, 

to be a simplified analysis, since that torque could also be computed 

indirectly as a byproduct of the next year's NASA-grant research (but that 

work, determining the actual wobble-induced flow in realistic oceans, is 


extremely difficult and could require more than a year's effort...)* Thus, 
the model considered here consisted of a single, vertical-sided continent in 
oceans of uniform density p and depth h overlying a spherical earth (radius " 
R); the flow was taken to be steady and unforced, and—slnce we're looking for 
non-dissipatlve coupling — the oceans were treated as inviscid. 

Prior to the analysis, fluid dynamics texts were consulted (Pedlosky, 

1979, chap. 1-3, 6, in which most of the concepts discussed here can be found; 
Tritton, 1977, chap. 1-15; portions of Lamb 1932, Batchelor 1968, Gill 1983, 
and a Phillips 1963 review article). The first step in the analysis was to 
scale the equations. Dimensional analysis quickly revealed that, because of 
the very slow long-period currents and global length scale characterizing pole 
tide-type flows, the situation being considered here corresponded almost 
exactly to geos trophy , where the dominant forces acting on the fluid are the 
Coriolis and pressure forces. I briefly investigated geostrophy on a 
sphere — a topic infrequently found in the literature due to the popularity of 
the beta-plane method. With no radial variations in any quantities and no 
radial flow, as appropriate for thin oceans, a geostrophic force balance 
allows only uniform zonal flow. If a continent is added to the oceans, then 
the presence of such a barrier creates a "Taylor band” of stagnant water in 
the same latitudes as the continent (this is the spherical equivalent of the 
Taylor column found in flat geometry). 

In order to treat realistic flow problems, then, either radial motion or 
departuc j from geostrophy, or both, must be allowed; for our situation the 
'force* most likely to break the geostrophic balance is advecticn, creating an 
inertial boundary layer mar the continent. The quasi-geostrophic continuity, 
momentum, and vorticity equations are, in vector form, 

V- u = o 
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where u * (u r> u tf , u^) * (0, u tf , u^) is the velocity and u> - V x u is its 
vorticity; and is earth's angular velocity. 


These equations were also scaled for thin oceans and then reduced 


accordingly. Functional analysis demonstrated that the only near-geos trophic 
solutions away from the continent were trivial ones, if the ocean depth h was 
not a function of longitude X* For the two-dimensional flow on a sphere, an 
obvious stream function T could be defined s - J, Uj ) 

and it was also shown that the radial vorticity could be written as 

= -£7 

this expression reveals two regimes for the flow, an "interior" nearly 
geostrophic flow (since <J r «v7 i Y) and a boundary flow (near the boundary ^h 
will be huge). For the model being considered, 11 is essentially constant 
within the oceans, yielding apparently an entirely geostrophic flow except for 
an infinite, infinitely thin inertial boundary layer. 

♦ 

These conclusions were unchanged even after the initial equations were 
re-scaled to allow explicitly for the possibility of a boundary layer of 
intermediate characteristic length. In a last attempt to uncover worthwhile 
'spherical* solutions, perturbation methods were tried. For the pole tide the 
Rossby number "c" is tiny, so all variables in the 'spherical' equations were 
expanded in powers of £ and then matched up. Unlike the flat earth case, the 
combination of continuity and momentum equations here is not "geostrophically 
degenerate" and instead led to 

U t ,50, U A * U x l\) > 7 * ?(\) 

where *? is the sea-level height; then, with a continental barrier present, the 
equations once again implied the existence of a Taylor band at the same 
latitudes as the continent. 

It therefore appeared necessary to resort to the well-known beta-plane 
approach if this brief research project were ’o be fruitful. On the 
beta-plane (centered at colatitude 9 m <9q), the scaled quasi-geos trophic 
equations were reduced for thin oceans. With L representing the horizontal 
length scale of the flow, then for the case where L/R is of order t the 
governing equation for the flow is (Pedlosky, §6.3) 

T( f>„ > + py) = O 

where J is the Jacobian operator, p Q is the lowest order non-hydrostatic 
pressure, ^ is a small parameter, j*> = tan0. , and y - (£(,-£ )R is the 


northward variable (x -JRainP* is eastward). 

Four situations were considered, in order to build gradually to the 
project goal (this Jacobian equation is highly non-linear! ). In all cases the 
topographic drag force D and drag torque T on the continent were calculated as 

JO * X + yM?. Uh* 

T - u*ii» 

integrated over the coastline/, which makes an angles with the eastward 
direction. 


I) uniform zonal flow u x - -0 from infinity, blocked by an infinite 
north-south continent (Pedlosky’s example of Inertial boundary layer). 
Results are 


where f 


0 
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2fLc.os9 t is the Coriolis parameter; and 

*0 T = 0 


(the drag is necessarily zero since the coastline is infinite). 


II) same as I but with a square, north-south/ jast-west continent whose 
sides are of length L. Results are 

po a L ^-^-7; expl-tp-r 

j5 = f * gt.vun'a-v. 

III) same as II but with flow at infinity similar to the pole tide. 

True pole tide currents on a spherical earth in global oceans — the "infinity" 
of the beta plane — would in beta-plane terms approximate to 

cotiO, r * 7 /* x 
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from Dlckoan (1984), where T * /jf aSilf '■**/§" j , ' T,*’ being the 
major spherical harmonic component of the pole tide in response to a wobble of 

amplitude M . However, the closest approximation to these u , u which is 
p x y 

consistent with the beta-plane equations was merely 

«* -ip/y/i .) 1 u, -r , 




with D ■ "l* . Even with this velocity at Infinity, barely more complicated 
than that In I and II, the shortcut method of solution employed by Pedlosky 
(see § 3.13) could not be used here. Postulating that the solution includt * 
exponential decay eastward, we eventually found 

p* ” -£-P.Ux - £ + &.*<■' . 

Conceivably the coefficients are given roughly by A ffl ■ 2p(L/2m «-) , . 

The net drag is 

/) = f-kUhL 1 £*[?//*-•] + y [2 ![&»,- }/l\ ^ 
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IV) same as III but with a circular continent c€ radius L. How the 

nonlinear Jacobian equation is complicated by highly non-linear boundary 

conditions. Although an analytic solution for the A . B as above is clearlv 
— m m J 

beyond the scope of this project, so that p„ remains unknown, It could be 
proved in this case that the net drag is independent of the Am, Bm. We found 

•4 

For flow with the characteristics of the pole tide, U * f 1.6x10 

8 

cm/sec; if the beta plane is centered at £: and L ^ 10 cm, then 

<>». [ 2.T> X + 0 . b y ) * I O ^ 

~f~ (l.S y ~ O.H x ) * l® , 
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With 7 such continents but slightly larger in size, L 5x10 cm (if the beta 
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plane approach is still valid for such L/R), the net torque would be ~ 3x10 
dyne-cm assuming the results for IV still apply; this is still several orders 
of magnitude smaller than required by Dickman (1983). 

It remains to be seen, when the goal of the next year’s project 
(Supplement 3) is achieved, whether the actual oceans— with irregula*’ 
continental coastlines, bottom topography, and unsteady turbulent flow— would 
exert couples on the solid earth <*10^ times greater.... 



My grant-related activities during the reporting period also included 
presenting a talk on the self-consistent dynamic pole tide at the Spring 1984 
AGU oeetings in Cincinnati (a copy of the abstract is attached); and the 
supervision of two graduate students (0. Steinberg, C. Ammon) . Steinberg has 
been working for over 1 1/2 years on the consequences of Dahlen's (1976; see 
also Merriaa 1973) self-gravitating, loading static pole tide model. 
Previously, the FORTRAN programs he had written to lmpleaent Dahlen's theory 
had yielded tide characteristics that did not agree with Dahlen's published 
values. During the reporting period, Steinberg discovered that the published 
values for the tide's effect on the Chandler wobble frequency laplldtly 
included the effect of the mantle's response to the pole tide load. This had 
not been clearly stated in the literature and explained why a tide supposedly 
enhanced (through self-gravitation and loading) compured to earlier 
non-loading static tide models would affect the wobble period less . 

Steinberg's programs were eble to separate the tide's own effect from that of 
the loaded mantle, and he showed that the loading, gravitating tide does 
indeed affect the wobble period more than the non-loading tide [however, the 
mantle response to loading reverses the net effect!]. It then became clear 
that the loading, gravitating tide i£ enhanced in magnitude, compared to 
non-loading pole tide models, at least according to the tide's spherical 
harmonic 2, -1 coefficient (which determines the tide's effect on wobble 
period). The next step — the final one in Steinberg's M.A. thesis 
research — will be to construct the loading, gravitating tide at various ports, 
and compare it with the actual pole tide; such a comparison (never before 
published) would allow departures of the actual tide from equilibrium to be 
accurately quantified. 

Ammon completed a year-long investigation of techniques for deconvolving 
Chandler wobble data. It was hoped that, by treating the Markowitz wobble as 
a second free wobble of the earth, the deconvolution filters of Smylie et al. 
(1973) could be modified to generate more accurate wobble excitation time 
series. Unfortunately the project got bogged down in reproducing the original 
Smylie et al. filters: perhaps that paper was a bit vague in providing de fc *ils 
of the filters; in any event, as constructed by us the filters appeared 
unstable, with undesirable side-effects. Without a much longer and Intensive 
effort, it is not likely that those filters could be "fine-tuned" enough to 


generate recognisable or reasonable excitation time aeries 


Following the Spring AGU meeting I reviewed once again the Supplement 1 
research upon which the talk was baaed. The review revealed several errors, 
primarily involving the normalisation used for the complex spherical harmonic 
functions. After correcting these errors it turned out that my theoretically 
computed self-consistent dynamic pole tide was, at least in global oceans, 
very similar to the idealised static pole tide; for example, the effect of the 
dynamic pole tide on the Chandler wobble period is likely to be within 1 day 
(l.e. 3Z) of the static tide's effect. I completed the manuscript on "the 
self-consistent dynamic pole tide In global oceans, M and submitted it to 
Geophysical Journal of the Royal Astronomical Society at the end of the 
reporting period; because of the greatly differing conclusions resulting from 
the revisions, I also circulated preprints to a number of colleagues* To save 
on photocopying and mailing costs, I will be mailing reprints to NASA. 
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The Self-Consistent Dynamic Pole Tide 


S.R, DICKMAN (Dept, of Geological Sciences, State 

University of New York, Binghamton, N.Y. 13901) 

It has been known since Newcomb (1892) that the pole 
tide, the oceanic response to the Chandler wobble, 
acts— through Its m?ss redistribution (products of 
inertia)— to lengthen the Chandler wobble period. The 
extent of lengthening, traditionally calculated by 
assuming a static response. Is ^ 1 month. 

Such estimates, however, fall to treat the pole tide 
self-conslstently. The tide's effect on wobble 
period, computed according to conservation of angular 
momentum (Llouvl lie equation), Is a function of the 
tide height; the latter Is found rigorously according 
to conservation of linear momentum (e.g. Laplace tidal 
eqns) and depends on t.he wobble period. The two must 
be determined jointly or the 1engthen1*q estimate will 
be Incorrect. 

At the Fall 1983 Meeting I reported on a first inves 
tigatlon of the self-consl stent pole tide. The oceans 
were taken as global, and linearized bottom friction 
of different strengths was considered. In all cases 
the resulting tide lengthened the wobble period by 
only % half the static amount. 

But a non-static tide must be supported by non-zero 
currents, and relative angular momentum associated 
with the currents will also affect the wobble period. 
Self-consistent theory Incorporating such effects has 
now been developed. Results to date indicate that the 
angular momentum contribution to wobble period length- 
ening Is surprisingly large: more than three times 
that of the tide's products of Inertia. 

Although oceanic non-globallty can be expected to 
modi fy these results, It appears that the self- 
consistent and dynamic nature of the pole tide is 
crucial to the tide's effect on wobble. 
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